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In this paper we summarize our STM studies of the density of electronic states in nearly optimally 
doped Bi 2 Sr 2 CaCu208+5in zero field. We report on the inhomogeneity of the gap structure, density 
of states modulations with four-lattice constant period, and coherence peak modulation. 



I. INTRODUCTION 

Tunneling spectroscopy has been an important tool 
in the study of high-temperature superconductors since 
their discovery. In the early days of high-T c a va- 
riety of gap sizes and structures were found and in- 
troduced much controversy into the subject. How- 
ever, recent measurements have been more consistent 
among groups, revealing a relatively coherent picture of 
the surface of high-Tc materials as viewed with STM 
CO S S El lJ • In this paper we review our work on 
the fine-scale structure of the electronic states at the sur- 
face of Bi?Sr?CaCuoOs^(BSCCO) as revealed by STM 
measurements 0, EJ ED] • Starting with an analysis 
of the shape of individual spectra, we further look at the 
spatial dependence of features such as the size of the gap, 
coherence peaks and local density of states (LDOS) mod- 
ulations. We perform measurements using a cryogenic 
STM on near optimally doped BSCCO (T c ~ 86K) grown 
by a floating-zone method. The samples are cleaved in 
an UHV of better than 1 x 10 -9 torr and then quickly 
lowered to the cryogenic section at a temperature of 6- 
8K. Most data were taken with a sample bias of -200mV 
and a set point current of -lOOpA At each point on the 
surface, dl/dV spectra was also taken. 



II. INDIVIDUAL SPECTRAL SHAPE 

A representative spectrum is shown in Fig. HJi. To an- 
alyze the spectrum we use a d-wave gap formula A(#) = 
Acos(20) with thermal broadening /c#T, and smearing 
T, averaged over all in-plane k-directions with a weight 
function g(0) appropriate for tunneling perpendicular to 
the Cu02 planes. 
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Where f(E — E') is a Fermi function and the Lorenzian 
broadening is sometimes replaced by a gaussian one with 
width T. In general a finite V tends to reduce the co- 
herence peaks. However, as can be seen from Fig. [Qd, 
spectra with small gaps tend to have unusually tall co- 
herence peaks. Thus, for the fit in Fig.QJi we use V — 0, 
obtaining A = 32mV and ksT = 0.7mV ± 0.3mV (con- 
sistent with the measurement temperature of 8 K). 
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FIG. 1: a) Typical spectrum with fits. Blue curve shows best d- 
wave fit, red shows d-wave fit with angle-dependent Fermi velocity, 
b) Series of spectra with coherence peaks for increasing gap sizes. 



The fit fails to capture elements of the data in several 
areas. First, there is always a significant particle- hole 
asymmetry in the data. While this probably comes from 
the non-flat shape of the normal state DOS, and could 
be added to the fit in the form of an asymmetric back- 
ground, the d-wave form, even with T = can still not 
fit the excessive height of the coherence peaks, nor the 
dip at approximately twice the gap size which is more 
pronounced on the negative bias side. The need for a 
g{0) that is different from unity was first proposed by 
Oda et at. [l3| to explain the more rounded shape of the 
gap at zero bias. The red curve in Fig. EJi shows a fit 
of this type, using g(6) which varies by a factor of two 
between the node and the antinode, yielding A = 31mV 
and ksT = 1.2mV ±0.3mV. This fit captures the shape 
of the subgap DOS, though there appears to be a slight 
excess of measured states at very low bias, perhaps due to 
zero bias anomalies. This fit accentuates the excess states 
at the coherence peaks as it fits the full amount of states 
removed from the gap region, but produces a coherence 
peak that is weaker than the one found experimentally. 
This effect is strong at small gaps and become weaker 
with larger gaps and lower coherence peaks. While there 
is no theoretical explanation for this result, this finding 
may point to a peculiar effect in which in the supercon- 
ducting state the coherence peak gains states from high 
energies as well as from the gap region. The fact that 
small gaps show larger coherence peaks may point that 
this effect is stronger in more overdoped samples as will 
be discussed in section IVTl For larger gap and lower co- 
herence peaks spectra a finite T is needed. A spectrum 
with a gap size A = 42 mV will typically yield T « 2 — 3 
mV if Eqn. ^ is used. While a precise determination of 
the numerical value of the gap is therefore impossible at 



2 



present, there is fairly good agreement between the gap 
values found by the fits and the location of the maximum 
in the dl/dV, so we will use this maximum G(A) as a 
phenomenological measure of A, throughout. For spec- 
tra that do not show coherence peaks we will use the edge 
where the conductance decreases below the background 
as a measure for the gap Q. 

The spectra shown in Fig. ^ are typical when coher- 
ence peaks are present. There are several observations we 
can make. First is that the gap can vary by as much as a 
factor of 2 within a scan area. Second, it is clear that the 
larger gaps have lower coherence peaks. Also, the energy 
of the dip feature seems to move with increasing gap size. 
This dip has been attributed to a strong coupling effect 
Q, regardless, it is apparently related to superconduc- 
tivity since its energetic location is scaling with the gap. 
Finally, the negative bias slope increases with gap size, 
while the bias asymmetry decreases. Presumably these 
latter changes reflect the evolution of the normal (non- 
superconducting) density of states. Also note that the 
asymmetry in the coherence peak heights changes with 
the gap size. The largest gaps have more weight in the 
positive bias coherence peak, while the smaller gaps have 
more weight in the negative bias peak. 



III. SPATIAL VARIATIONS OF THE 
SUPERCONDUCTING PROPERTIES 

Fig. [2ji shows a typical map of the size of the gap. 
Patches about 30 A across of varying gap size are clearly 
visible. While the magnitudes of the largest gaps may 
vary between samples (mostly dependent of fraction of 
gaps with no coherence peaks), the magnitude of the 
smallest gaps observed is always ~ 30 mV. The smallest 
scale features reflect some variation with atomic resolu- 
tion, and the partial near vertical lines show that there 
is some correlation between superstructure and the gap. 

Spatial variations of the gap on the surface of BSCCO, 
similar to Fig. [2Jl were reported by several groups 
0O Ell However, the origin of the inhomogeneities is 
still not clear. Martin and Balatsky proposed a phe- 
nomenological model where disorder in the dopant sites 
will lead to variations in the carrier density in the plane 
due to the unusually large screening length of doped Mott 
insulators. The local doping is therefore the average den- 
sity of dopant atoms within a superconducting coher- 
ence length lead to patches with differing superconduct- 
ing properties, particularly the gap size. Wang et al. [l8| . 
used a t-J model to show that the poor screening results 
in an effective, larger than the actual distance, separation 
between the Cu02 planes. Both calculations showed that 
simultaneous measurements of the spatial extent of the 
doping variations and the accompanying gap size varia- 
tions provides a check on whether the doping inhomo- 
geneities can account for the gap inhomogeneities. 

Fig. [3| shows a histogram of the gap size over the area 
shown in Fig. [2t with the best gaussian fit to this his- 
togram. The width of the gap distribution is about 7.5 




FIG. 2: Comparison of a) gap size and b) coherence peak height 
over a 140A x 140A area. Arrows denote the Cu-O bond direction. 



mV. The model of Martin and Balatsky, using the ex- 
perimental value of the coherence length, ~15A, and a 
gap that is approximately linear in doping with a coeffi- 
cient of 0.3V/carrier, yields a standard deviation of the 
gap size of ~30mV . This is much larger than we found 
experimentally. The model of Wang et al. uses a dif- 
ferent gap versus doping dependence but yields similar, 
though slightly smaller, discrepancies. The above analy- 
sis may point that other models that are more qualitative 
and based on electronic phase separation may be the true 
explanation for the gap inhomogeneities phenomenon as 
was previously discussed by Howald et al. 0. 

The spatial variation of the coherence peak height 
differs somewhat from the spatial variation of the gap. 
Fig. [213 shows the variation of the coherence peak height 
in the same area of Fig. While there is a definite 
correlation between gap size and coherence peak height, 
there are several differences. First, the peak height varies 
considerably on the atomic scale, while the peak position 
does not. Second, the peak position more clearly exhibits 
the granular structure: except for the atomic scale vari- 
ations, the variation in the peak height is smoother. Fi- 
nally, it is clear from Fig. that the peak-height shows 
an ordered structure, especially if one considers the lower- 
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FIG. 3: a) Histogram of gap size as dervived from Fig. |2Jl. b) 
Coherence peak conductance vs. gap size derived from Fig. |2] 



right corner of the figure. We discuss this phenomenon 
in section IVll 



IV. PERIODIC STRUCTURES IN THE LDOS: 
GENERAL CONSIDERATIONS 

The discovery [HI [2(J of stripe order in 
Lao-^Sr^NiO^^ and soon after in Lai.6- £C Ndo.4Sr, E Cu04 
|22 | added considerable credibility to the suggestion that 
charge ordered states form an important bridge between 
the Mott insulator, and the more metallic state at heavy 
doping. STM however is a static probe and thus can- 
not detect any structure associated with fluctuating or- 
der unless something pins it. Indeed the inhomogeneities 
discussed above and other point defects are a natural 
source for pinning. As a result charge order will be visi- 
ble to STM in the form of LDOS modulations. However, 
defects can also create other effects that need to be un- 
derstood before a static charge modulation explanation 
is invoked. 

Defects in simple metals will also cause LDOS mod- 
ulations. These "Friedel oscillations" |23| are related to 
Fermi-surface-derived non-analyticities in the suscepti- 
bility, x(k). A generalized form of these oscillations can 
occur in more diverse systems in which the relevant struc- 
ture in x(k) is not directly related to any feature of a 
Fermi surface. In particular, if the system is proximate to 
a charge modulated state, such as a stripe state, the val- 
ues of k = q at which x nas maxima will reflect the pat- 
tern of spatial symmetry breaking of the ordered state, 
but x(k) will respect the full point-group symmetry of 
the crystal (unless the liquid state is a nematic fl2|). 
So, the generalized Friedel oscillations around a point 
impurity in a stripe-liquid phase will inevitably form a 
checkerboard pattern |l2|, |24| . 

There is another form of spatial modulation of the den- 
sity of states, one with a period which disperses as a func- 
tion of the probe energy. This latter effect, which was 
first demonstrated by Crommie et al. 25], is produced 
by the elastic scattering of quasiparticles of given energy 



off an impurity. The resulting interference between scat- 
tered waves leads to variations of the local density of 
states at wave vectors q = k — k', where k and k' are 
the wave- vectors of states with energy E^ = = ek', 
as determined by the band structure, Generalized 
versions of these oscillations can occur even when there 
are no well defined quasiparticles, so long as there are 
some elementary excitations of the system with a well- 
defined dispersion relation. However, quasiparticle scat- 
tering interference will take place at a certain energy and 
wave- vector only if there are available states. Concen- 
trating on q^-o = |(27r/ao), It is easy to see that the 
quasiparticle interferen ce p icture cannot produce such 
peaks at low energies [13, |2£|. Assuming quasiparti- 
cles with k-dependent energy: = yj A^. + where 
Ak = ^-\cos(k x a§) — cos(k y ao)] is the d-wave super- 
conducting gap. For the above q vector we can take 
In that case k = (— 7r/4ao, k y ), and 
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q = (2^/4a ,0) 
k' = (7r/4a ,fe 
ARPES data h 
A ~ 30 meV |2 




If we extract k y directly from the 
I23l, we find k y ~ 0.6(7r/ao). For 
U, which is also the minimum 
gap found on this type of samples using STM, this 
gives an estimate of the lowest energy for quasiparticle 
scattering at this wavevector of approximately Aq/2 ~ 15 
meV. This energy cutoff could differ a little, depending 
on the exact details of the band-structure, but this anal- 
ysis certainly excludes the energies around zero bias. 

However, the above sources for LDOS modulations 
should not appear exclusively. In fact, charge order and 
quasiparticle scattering are likely to coexist in a system 
with well defined quasiparticles. 0, |^ IHH The rela- 
tionship among these effects will be discussed below. 



V. PERIODIC STRUCTURES IN THE LDOS: 
EXPERIMENTAL RESULT 

The essence of the above discussion is that in STM 
studies of cuprates we would expect stripe or checker- 
board correlations to make an appearance as generalized 
Friedel oscillations, while quasiparticle-like interference is 
a distinct phenomenon that could also be present. The 
observation of a checkerboard pattern with a ~ 4ao pe- 
riod about vortex cores in BSCCO |34( has been a pos- 
sible evidence for pinned charge stripes. However, it re- 
mained of great importance to find evidence of charge 
modulation in BSCCO with no applied field. 

Indeed, Howald et a/.|9[ shortly afterwards reported 
this same effect in zero field on similarly doped BSCCO 
crystals fabricated without intentional substitution of 
impurities. The observed modulation with ordering 
wave vector q^-o ~ [0.25 ± 0.03] (27r/ao) was found at 
all energies, exhibiting features characteristic of a two- 
dimensional system of line objects. Moreover, Howald et 
al showed that the LDOS modulation manifests itself, 
for both positive and negative bias, as a shift of states 
from above to below the superconducting gap. The fact 
that a single energy scale (i.e. the gap) appears for both 
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superconductivity and these modulations suggests that 
these two effects are closely related. The summary of 
the results of Howald et al. is shown in Fig. 0] empha- 
sizing the fact that peaks in the Fourier transform at 
q^-o ~ [0.25 ± 0.03](27r/ao) are found even at low en- 
ergies where quasiparticle scattering interference should 
not produce a signal. Fig. 0] was reproduced theoretically 
by Podolsky et al. |3^ . 
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FIG. 4: Fourier transform at q=(27r/4ao)(0,±l), the location of 
the peaks as a function of sample bias. The red and blue traces 
correspond to the real and imaginary parts, respectively. Right: 
Two Fourier transforms maps at low energies showing peaks at 
qiT -o ~ [0.25 ± 0.03] (27r/a ) and qo-^ ~ [0.25 ± 0.03] (2ir/a ). Cir- 
cles denote the main contribution to the peaks. 

Subsequent studies at zero field [3^, |3(j measured the 
dispersion of the strongest Fourier peak along the tt — 
(i.e. Cu-O) direction. They asserted that it was con- 
sistent with what is expected from quasiparticle scatter- 
ing interference. [26| In general, their data showed good 
agreement with photoemission results (i.e. band struc- 
ture results [37}) at large bias, but unlike photoemission 
results, did not continue to disperse below ~ 15 mV. 

To resolve the discrepancy Kivelson et al. ^3 P ro ~ 
posed several procedures which allowed Howald et al. to 
separate the two effects. First, integration of the LDOS 
over a wide range of energies reduces the influence of any 
random or dispersing features such as quasiparticle scat- 
tering interference, while at the same time it enhances 
features that do not disperse. Fig. [5| shows the result of 
such a procedure. 

A complementary approach to separate pure charge 
modulation effects is to look for an interaction of the 
charge modulation with the superconducting order pa- 
rameter. We claim that the periodic structure observed 
in G(A) in Fig. 03 is exactly this effect. We discuss this 
procedure in the next section. 



VI. MODULATION IN G(A) 

Coming back to Fig. we compare the coherence peak 
heights and the gap size. The data for this figure was 
taken with a setpoint voltage of 65 mV. The pronounced 
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FIG. 5: Line scans as a function of k x along the (0,0) to (tt, 0) 
direction, and as a function of energy (color scale). Top panel shows 
the LDOS (dl/dV), and bottom panel the integrated LDOS (I) up 
to the given energy. 



enhancement of the modulation signal with this setpoint 
voltage led us to a novel procedure to normalize (i.e. di- 
vide) the individual spectra by the current at +65 mV. 
A detailed description of this procedure is given in Fang 
et al. [ll|. Fig. |Hl shows a representative Fourier trans- 
form of low bias together with a Fourier transform of the 
peak in the dl/dV (i.e. G(A)) for a 160A x 160A area 
sample |9|- A clear correspondence of the peaks in the 
— 7r direction is found when comparing these two maps. 
This is a striking result since A is in the range of 30 - 
60 mV which according to band structure is the strongly 
dispersive region of quasiparticles scattering interference. 

Fig. El also shows a line scan in the — ir direction 
together with G(A) for the same sample. For all our 
samples, our Fourier analysis from low energies up to the 
smallest gap sizes (where the noise from inhomogeneities 
overwhelms our signal) supports the picture presented 
earlier of a non- or weakly dispersive feature in the re- 
gion q = 0.22(27r/ao) — 0.25(27r/ao) in addition to disper- 
sive features at a lower k- vector. We find that the large 
amplitude of the lower /c-vector features swamp out the 
non-dispersive feature at higher energies, but become rel- 
atively weak at lower energies. The additional line scan 
for the Fourier transform of coherence peak heights can 
be seen as a way to remove the effects of gap size inho- 
mogeneities to reveal that a structure at q w 0.25(27r/ao) 
still exists at higher energies. It makes a similar point as 
the spatial maps of the coherence peaks, namely, that by 
selectively sampling from the (higher) energies related to 
superconductivity, the low energy features reappear. 

By comparing Figs.[2t andj^t), one can see that the am- 
plitude of the coherence peak DOS modulations is larger 
in the regions of large gap. In contrast, regions of small 
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FIG. 6: a) FFT of LDOS at 10 mV and of LDOS taken at the 
coherence peak maximum with 65 mV normalization (see text.) b) 
Dispersion relation of the charge modulation periodicity. Black line 
is normalized coherence peak maxima. 



gap show modulations of reduced amplitude. Since there 
are only a few modulation crests and troughs within a 
particular region of large or small gap, this effect is dif- 
ficult to quantify, although it can most easily be seen by 
following the regions of largest gap. We note that the 
regions of large gap with low coherence peaks generally 
resemble slightly underdoped samples Q. On the other 
hand, the modulation is suppressed and the coherence 
peak heights are more uniform in regions of small gap 
and tall coherence peaks. Gaps in these regions are more 
similar to gaps found in overdoped samples j^; this is 



consistent with the notion that beyond optimal doping a 
more homogeneous charge density exists closer to a Fermi 
liquid state. 

Such an observation does not necessarily point out 
a competition between charge-density modulation and 
superconductivity, but rather reinforces the idea that 
the two effects coexist at and below optimal dop- 
ing. One possible interpretation is that the fluctuating 
stripe/checkerboard phase exists in all the regions be- 
low optimal doping, and as one moves further into over- 
doping (i.e. into regions of small gap) the modulations 
become diminished. Our observations therefore comple- 
ment those of Vershinin et at. [^U who found similar pat- 
terns in the pseudogap regime of slightly underdoped 
Bi 2 Sr 2 CaCu208+5. As noted by Kivelson et al. 12], 
the effect of quasiparticle scattering interference should 
disappear at temperatures above T c , revealing the under- 
lying order. For our measurements at low temperature, 
in regions of very large gap with weak coherence peaks 
(similar to the pseudogap), charge ordering is indeed visi- 
ble. The two results therefore suggest that in the absence 
or suppression of superconductivity, charge-ordering may 
be the preferred phase. 



VII. CONCLUSIONS 

In this paper we discussed the evolution of structures 
in the LDOS on the microscopic scale. We first analyzed 
the individual spectra pointing to anomalies in their 
shape and discussing the determination of the size of 
the superconducting gap. Maps of the gap lead to an 
inhomogeneous pattern that cannot be explained by 
simple disorder in doping and thus may point to intrinsic 
electronic phase separation. The inhomogeneities also 
reveal ordered checkerboard patterns in the LDOS with 
a period close to four lattice constants. Studying the 
interplay between the two effects we suggest that both 
have a common origin. 
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